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Abstract 

A method is introduced to analyze the degeneracy properties of the band structure of a photonic 
crystal making use of the super cells. The band structure associated with a super cell of a photonic 
crystal has degeneracies at the edge of the Brillouin zone if the photonic crystal has some kind of 
point group symmetry. Both E-polarization and H-polarization cases have the same degeneracies 
for a 2-dimensional (2D) photonic crystal. Two theorems are given and proved. These degeneracies 
can be lifted to create photonic band gaps by changing the transform matrix between the super 
cell and the smallest unit cell which is caused by changing the translation group symmetry of the 
photonic crystal (the point group symmetry of the photonic crystal may keep unchanged). The 
existence of the photonic band gaps for many known 2D photonic crystals is explained through the 
degeneracy analysis. Some structures with large band gaps are also found under the guidance of 
this degeneracy analysis. 
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I. INTRODUCTION 



Photonic crystals, which are periodic arrangement of dielectric or metallic materials, have 
attracted a wide attention recently in both physics and engineering communities in view of 
their unique ability to control light propagation 0, |^, ^ . Many potential applications of 
photonic crystals rely on their photonic band gaps (PBGs). It is thus of great interest to 
design photonic crystals with an absolute band gap (as large as possible), particularly, for a 
given dielectric material. 

Two-dimensional (2D) photonic crystals have attracted special attention since they are 
easier to fabricate. Many 2D photonic crystals with large absolute band gaps have been 
found 1^, 1^, 0. Based on the difference of the filling factors between the dielectric band 
and the air band (related to the distribution of the displacement field), the rule of thumb 
sometimes can be used to explain the band gaps particularly at low frequencies ||. Due 
to the complication of the differential operators in electrodynamics (different field compo- 
nents are coupled to each other even if the permittivity e{r) is separable), it is difficult to 
obtain analytical (even approximate) solutions for the distribution of the displacement field 
(particularly at high frequencies). Therefore, many photonic crystals with large absolute 
band gaps can not be explained or found by the rule of thumb [|10|, [11], |15|. 

Degeneracy lifting is another explanation for absolute band gaps and even a method 
to create band gaps |]T2|, |T3|, [1^, |T^. The degeneracy can be lifted by e.g. using 

hexagonal photonic structures [Q, using anisotropic materials |p!3| , |T^, breaking the space 
group symmetry []T3|, [TB[ and changing the dielectric distribution without breaking the space 



group symmetry [|T^. Both accidental and normal degeneracies can exist in a photonic band 



structure (see e.g. [T^], this is different as compared to an electronic system). To investigate 
the degeneracy properties of 2D photonic crystals, the i?-polarization and i/-polarization 
are usually considered separately as suggested by |]12|, ITSf. It is very complicated to predict 
where the degeneracy appears and how to break the degeneracy. 

In some cases, we don't have to rely on such an analysis. In the present paper, we 
introduce a method to create degeneracies first and then break them to create band gaps 
by studying the band structure associated with a super cell (instead of the unit cell as 
considered in the literature mentioned before by others). In the band structure associated 
with a super cell, we can analyze how degeneracies are formed and how to break them to 
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create band gaps. 

The point group symmetry of a photonic crystal is defined with respect to the point 
with the highest symmetry. For example, the point group symmetry does not changed by 
adding columns at the corners of the unit cell for the 2D photonic crystals considered in 
(they belong to the same point group symmetry C^^). We notice that the translation 
group symmetry does not change either. Thus the space group symmetry of the photonic 
crystal does not change at all although the smallest unit cell must include two rods after the 
additional rods are added. It may be hard to understand the degeneracy breaking for the 
iJ-polarization at point M of the second and third bands as shown in |jl5| without careful 
analysis of the electromagnetic field distribution. However, if we study the band structure 
associated with a super cell, the lifting of the degeneracy and the creation of PBGs of such 
photonic crystals can be understood with some tricks even when the space group symmetry 
of the photonic crystal does not change. For the example mentioned above, the photonic 
crystal with the additional columns at the corners can be treated as the result of changing the 
translation group symmetry (keeping the point group symmetry unchanged) from another 
photonic crystal with the additional columns having the save size with the original column 
(as shown in Fig. 3(a) with the square column case), which also belongs to C4i, point group 
symmetry ). Another view is presented to understand the idea of additional columns. Not 
surprisingly, many known photonic crystal structure such as a chessboard structure |^, a 



square lattice of square rods [|n|] and even a triangular air hole structure can be somehow 
understood from this point of view. (cf. the numerical example associated with Figs. 3-6 
below). By using such a degeneracy analysis associated with a super cell, some structures 
with large band gaps are also found in the present paper. 



II. THEOREMS FOR DEGENERACIES IN THE BAND STRUCTURE ASSOCI- 
ATED WITH A SUPER CELL 

A unit cell we consider here refers to the smallest periodic region in a photonic crystal. 
If the periodic region includes more than one unit cells, e.g. 2 unit cells, it is called a super 
cell. First we want to study the relation between the band structure associated with the 
super cell and the original band structure (associated with the unit cell). 

In general, we consider a three-dimensional (3D) photonic crystal with primitive lat- 
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tice vectors ai, a2 and as. The associated primitive reciprocal vectors bi, b2 and bs are 
determined by 

3 

b. = 27:^^ -, (1) 

ai • (a2 X ag) 

where tijk is the 3D Levi-Civita completely antisymmetric symbol. The complete set of 
the reciprocal lattice vectors is written as {G|G = /ibi + l2h-2 + /sba}, where (/i, /2, are 
integers. We denote the first Brillouin zone formed by these reciprocal lattice vectors {G} 
as zone A. 

The primitive lattice vectors for a super cell are the linear combinations (with integer 
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coefficients) of the primitive lattice vectors for the unit cell, i.e., a'j = E ^ij^ji hj — 1; 2, 3, 
where Nij are integers. The corresponding primitive reciprocal vectors for the super cell are 

3 

determined by b^ = '^^ ^'''2'^.{a! xa!) ■ '^^^ integers N^j form a 3 x 3 transform matrix with a 
positive determinant det(A'") = M > 0. 
Since 
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b. • a; = 2n^^^ — • E 7V,,a, 

ai • [a.2 X aaj 



3 

z 

1=1 



(2) 



it follows from a'j • h'j — 2'KSij that 



b. = E%b;. = E^5b;, (3) 
j=i 3=1 

where the superscript T denotes the matrix transposition. The set of the reciprocal lat- 

3 3 
tice vectors associated with the super cell is {G'|G' = E ^j^'j}- Since G = E f^ihi = 

j=i i=i 

3 

J2 i^iNljh'j, one sees that {G} is a subset of {G'}. Note that the elements of {G} and 
{G'} are the integer grid points (they do not fill any continuous space) formed by the 
corresponding reciprocal lattice vectors. We denote the first Brillouin zone formed by the 
reciprocal vectors {G'} as zone B. 

Lemma. There exist a subset {G} of {G'}, which satisfies: 

(i) {G} C {G'} and {G} n {G} = 0. 

(ii) There are M elements in the set {G} (M is the determinant of the matrix N) and 
the difference of any two of them does not belong to {G}, i.e., (Gi — G2 ^ {G}). 



(iii) Any G' G {G'} can be expressed as 

G' = G + G, (4) 

where G G {G},G G {G}. 

The proof and a way to find the set {G} are given in the appendix. 

If we define the addition of vectors as the multiphcation in a group theory, we can take 
{G'} as a group and {G} as a subgroup. Then the vector is the unit element of the group. 
From the group theory, one knows that {G'} is the union of all the cosets of the set {G}. 
The subset {G} is used to give the cosets. 

With these reciprocal vectors, each eigenstate of the electromagnetic field component ifk 
(with the wave vector k in the first Brillouin zone) in the photonic crystal can be expressed 



in terms of the following Bloch series [10 



iJk(r) = e^^-^^i/Ge^^^ (5) 



G 

The field component ifk satisfies the following equation, 

e/7k = (6) 

where the operator can be easily derived from Maxwell's equations, and c is the speed of 
light. 

For any wave vector k in the k space of a photonic crystal, one can find a wave vector in 
the first Brillouin zone which has the same eigenstate. The difference of the two wave vectors 
should be a reciprocal vector. Therefore, for any wave vector k, there exist a G G {G} so 
that 

ki = k - G (7) 

is in zone A (associated with the unit cell) and a G' G {G'} so that 

k2 = k - G' (8) 

is in zone B (associated to the super cell). We call ki (in zone A) the counterpoint of k2 (in 
zone B) for the same photonic crystal. They denote the same eigenstate in the reciprocal 
vector spaces associated with the unit cell and the super cell, respectively. 
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For a fixed k2 G B, we define the set {Ki|Ki = k2 + G, for all G e {G}}k2- Clearly, 
there are M elements in {Ki}. Since not all of these M elements are in zone A, we can make 
each of them inside zone A by subtracting an appropriate reciprocal vector Gki G {G}. 
Thus, we define a set |ki|ki = Ki = Ki — Gki in A for all Ki G {Ki}k2}k2- Obviously 
{ki}k2 contains M points inside zone A. 

Theorem 1. 

(i) The M elements in {ki}k2 are the counterpoints of k2. They are M different points 
in zone A. 

(ii) All the M eigenstates with the M wave vectors in {ki}k2 (associated with the unit 
cell) correspond to M eigenstates with one wave vector k2 in zone B (associated with the 
super cell). 

(iii) Each band in the band structure associated with the unit cell will split into M bands 
in the band structure associated with the super cell. 

Proof. For any wave vector k in the k space, ki and k2 are the counterpoints in zone A 
and zone B, respectively. From Eqs. (4), (7) and (8), one has 

k = k2 + G' = ki + G, (9) 
ki = k2 + G' - G 
= k2 + G + Gi - G 
= k2 + G + G2 

= Ki + G2. (10) 

Since ki is in zone A, it follows from the definition that ki G {k-^jkj (here Gki = — G2). 
Therefore, for any wave vector which has a counterpoint k2 in zone B, its counterpoint 
in zone A must belong to {kj^}k2- On the other hand, all the elements in {ki}k2 for all 
possible G' (correspond to all possible G G {G}) in Eq. (10) are all the counterpoints of 
k2. Therefore, the elements in the set {ki}k2 are exactly all the counterpoints of k2. 

Consider two different Gi, G2 G {G}. Correspondingly, we have Ki = k2 + Gi and 
K2 = k2 + G2. From the definition we have Ki — K2 = Ki — Gki — (K2 — Gkz) = 
Gi - G2 - (Gki - Gk2)- Since Gi - G2 ^ {G} and Gk^ - Gkj e {G}, we know that 
Gi — G2 and Gki — Gk2 are different, i.e., Gi — G2 — (Gki — GK2) 7^ 0, which immediately 
gives Ki — K2 7^ 0. This proves that Ki and K2 are two different points in zone A. Therefore, 
the elements in {ki}k2 are M different points. 
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Let Hki be the eigenstate for a wave vector in the set {ki}k2 associated with the unit 
cell. Prom Eqs. (5) and (10), one has 

G 

^ gi(k2+G)+G2-r J2 ifce*^ "" 
G 

^ gjk2T ^ ^^gi(G+G+G2) r 
G 

= e^^-^Y.HG'e''''-^HUr), (11) 

G' 

where H^^^ is the same eigenstate (with the same field distribution) but for the wave vector 
k2 in zone B (associated with the super cell). Therefore, all the M eigenstates with the M 
wave vectors in {ki}k2 (associated with the unit cell) can be represented by M eigenstates 
with one wave vector k2 in zone A (associated with the super cell). The M points of ki on 
any band in the band structure associated with the unit cell will be on M bands for one k2 
value in the band structure associated with the super cell. Generally speaking, one band in 
the band structure associated with the unit cell will split into M bands (which may overlap 
partially and form degenerated eigenstates) in the band structure associated with the super 
cell. The theorem is thus proved. 

Theorem 2. If a photonic crystal has some kind of point group symmetry, the eigenstates 
at the edge of the first Brillouin zone B will be degenerated in the band structure associated 
with the super cell. The degree of the degeneracy depends on both the determinent M of 
the transform matrix and the point group symmetry of the photonic crystal. 

Proof. For a wave vector k2 at the edge of zone B, besides that itself ki = k2 + 
(corresponding to Gi = 0) is one of its counterpoints at the edge of the zone B, it may 
have other counterpoint k^ = Ki — Gki (with Ki = k2 + G2) located at somewhere else 
at the edge of zone B. For nonzero G only those points at the edge of zone B may have 
counterpoints still at the edge of zone B and the counterpoints for those points inside zone 
B will be outside zone B (but still inside A according to the definition for counterpoints; 
note that zone B is inside zone A). Wave vectors ki and h[ correspond to the same wave 
vector k2 in zone B associated with the super cell. Sometimes there exists a symmetric 
operation a (which can be represented by a matrix for coordinate transformation; then 
one has — a^) and the associated operator T{a) (with T{a)f{r) — f{a~^r)) for the 
photonic crystal such that a\s.[ — ki and T{a)Q{r) = Q{a~^r)T{a) = Q{r)T{a). Assume 
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that i?ki and i?k'i are the eigenstates for these two wave vectors, i.e., ©iJki 
QH^', = ^H^',. Since 

T{a)H^>,{v) = //k'i(«-'r) = //,k'i(r) 

we have 

Gi/ki = er(Q;)i/k'i = r(Q;)ei/k'i 

= ^T(a)i/ka = ^i^ki, (13) 

we have a;ki = <^k'i- Therefore, ifki and ifk'i are two different eigenstates (for different 
wave vectors ki and k'l) with the same eigenvalue. In the band structure associated with 
the super cell, these two eigenstates are located at two bands but have the same wave vector 
k2 and the same eigenvalue. Thus, they are degenerated states. Since we do not assume any 
specific form for in the above proof, the theorem is valid in any dimensional space (and 
for any polarization in the 2D case). 

In the next section, we will illustrate these degeneracy theorems with some numerical 
examples and use the degeneracy analysis to explain PBGs for some known 2D photonic 
crystals and create large band gaps by breaking the symmetric properties of the photonic 
crystal. 

III. NUMERICAL RESULTS 

First we give a one-dimensional (ID) example. Fig. 1 is a ID photonic crystal consisting 
of alternating layers of materials with two different dielectric constants (si = 13 and £2 = 1)- 
We can select a periodic region (a super cell) to include 2 unit cells as shown in Fig. 1(a). 
The band structure associated with the unit cell and the band structure associated with the 
super cell (with N = 2) arc given in the same figure (see Fig. 2(a)), where the frequency 
and the wave vector are normalized with the same constant a = 1 in order to make them 
comparable. For this case, we have {G} = {Gi, G2} — {0, 0.5(27r/a)} and M — 2. From Fig. 
2(a) one sees that the eigenvalues (associated with the original unit cell) for the wave vectors 
outside the first Brillouin zone B (associated with the super cell) have their counterpoints in 
zone B in the band structure associated with the super cell. As expected, each band (solid 
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= ^H^, and 



(12) 



FIG. 1: A ID photonic crystal consisting of alternating layers of two different materials, (a) The 
super cell (including 2 unit ceUs) of the photonic crystal, (b) The symmetry of the super cell 

is broken (by changing the widths of the two dielectric layers while keeping the position of both 
unchanged) to form a unit cell for a new photonic crystal. 

line) associated with the unit cell corresponds to two bands (dashed lines) associated with 
the super cell. Since the center point of the super cell in Fig. 1(a) is mirror symmetric, 
one has a = —1. At the edge of zone B, the two wave vectors k[ = 0.25(27r/a) and 

ki = — 0.25(27r/a) — ak[ correspond to the same wave vector ^2 = — 0.25(27r/a) (note that 
^2 + ^1 = ki,k2 + G2 = k'-^. Thus, these two eigenstates are degenerated in the band 
structure associated with the super cell. Each eigenstate at A: = ±0.25(27r/a) is formed by 
two degenerated statcxs in the band strnctnre associated with the super cell. If one breaks 
the point group symmetry with respect to the center point of the super cell of the photonic 
crystal by changing the size of the inclusion media, one obtains a new photonic crystal as 
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FIG. 2: The corresponding band structures of the ID photonic crystals with si = 13, £2 = 1- (a) 
The sohd hnes are for the band structure associated with the unit cell and the dashed lines give 
the band structure associated with the super cell. Here we choose b = 2a and d = 0.5a for Fig. 
1(a). (b) The solid lines give the band structure for a new photonic crystal with a unit cell shown 
Fig. 1(b) (here we choose di = 0.3a and ^2 = 0.7a). The dashed lines are for the band structure 
for Fig. 1(a) before the symmetry of the super cell is broken. 

shown in Fig. 1(b). Since the resulted photonic crystal is still mirror symmetric respect to 
the center point of the inclusion media, the symmetry breaking with respect to the center 
point of the super cell does not change the point group symmetry of the photonic crystal 
as a whole. However because of the symmetry breaking, the translation group symmetry 
changes which leads to a larger unit cell. The corresponding transform matrix between the 
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FIG. 3: The case for the square lattice of square dielectric rods, (a) The super cell including 4 
unit cells, (b) The symmetry of the super cell is broken as two square rods increase in size and 
the other two rocls decrease in size . The symmetry- broken super cell is marked by the thick solid 
lines, (c) The first Brillouin zone A (marked by the dashed lines) associated with the original unit 
cell and the first Brillouin zone B (marked by the sohd lines) associated with the super cell of Fig. 
3(a). The first Brillouin zone C associated with the new unit cell (marked by the dashed lines of 
Fig. 3(b)) of the new photonic crystal of Fig. 3(b) is marked by the dotted lines. F, X, J, M are 
the symmetry points. 

super cell and the new unit cell changes from TV = 2 to A'^ = 1. The band structure is shown 
by the solid lines in Fig. 2(b), where one sees that the degeneracy disappears (since M = 1 
for this new photonic (crystal and consequently there is only one counterpoint for each wave 
vector in the Brillouin zone) and more band gaps appear. 

For the 2D case, if the dielectric inclusions have rectangle shapes, we can employ the 
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FIG. 4: The band structures associated with Fig. 3 for / = 0.47, e = 8.9. The sohd lines denote 
the H-polarization case and the dashed Unes are for the E-polarization case. They are calculated 
by the plane wave expansion method (with the inverse rule) with 289 plane waves for (a) (3 = 0, 
i.e., a simple square lattice of square rods; (b) /? = 1, i.e., the case of chess board. The absolute 
band gap Au/uc = 0.07 appears at Uc = 0.605^^7rc/a). (c)The bandgap map for < /3 < 1. 



plane wave expansion method with the inverse rule pTl to calculate the band structure. It 



is shown in ||2T| , that this method with 225 plane waves can give more accurate results 
than a conventional plane wave expansion method with even 1681 plane waves. In our 
calculations, we use this method with 289 plane waves and the error in the band structure 
is less than 0.5%. 

^2 



cell. The point group symmetry is C4„. We use the two mirror symmetries a 




As an example, we choose 2x2 unit cells as the super cell. Then we have A^^ 

"^0 2 

which corresponds to {G} = {0, b'l, b'2, b'l + b'2} (see the Appendix for a derivation for 
a general case). Fig. 3(a) shows a simple square lattice of square dielectric rods. Alumina 
is chosen as the dielectric medium and thus e = 8.9. The filling factor is set to / = 0.47. 
Fig. 4(a) gives the band structure associated with the super cell. Each band (associated 
with the unit cell) have split into 4 bands in the band structure associated with the super 

-1 
1 

and a2 = I I to analyze the degeneracy here. The eigenstates with wave vectors at 

\o -1) 

the 4 corners (e.g. points X and M) of the Brillouin zone B have 4-fold degeneracy and the 
eigenstates with wave vectors on two opposite edges are of 2-fold degeneracy under these 
two symmetry operations. Thus, in the band structure associated with the super cell one 
can see that there are 4-fold degeneracies at X and M points There are also 2-fold accidental 
degeneracies for the 4 split bands (and thus one sees only 3 bands) in region T — X. The 
-E-polarization and the if-polarization have similar behaviors of degeneracy. 

In the super cell, both the size and the position of the inclusions can be changed to 
break the point symmetry with respect to the center point of the super cell. Since the band 
structure is more sensitive to the inclusion size the size of the inclusions will influence 
significantly the band structure. In our first example, the symmetry is broken as two square 
rods increase in size and the other two rods decrease in size in order to keep the filling factor 
/ = 0.47 unchanged. In the resulted photonic crystal shown in Fig. 3(b), the squares rotate 
a 45° angle to form a chessboard structure in the new unit cell denoted by thickened lines 
in Fig. 3(b) after the symmetry is broken. The ratio of the side lengths between the smaller 
rods to the larger rods is 1 — /9 with < (3 < 1. When < /9 < 1, it just the case with the 
smaller square rods being added at the corner of the simple unit cell. When (3 = 1, the side 
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length of the smaller rods is and only two larger rods exist in the super cell. The structure 
is exactly the chessboard structure reported in [pO| . 

For < f3 < 1 shown in Fig. 3(b), since the symmetry with respect to the center point of 
each square (which has the highest symmetry) is still C4„, the point group symmetry of the 
photonic crystal is still (74^,. Similar with ID case, the symmetry breaking with respect to 
the center point of the super cell changes the translation symmetry of the photonic crystal. 
The corresponding transform matrix between the super cell and the unit cell changes from 

iV^ =1 I to A^^ = I I since the unit cell changes to a larger one. Therefore, 

[02) [o-i) 

the degeneracies of the band structure associated with the super cell will also change. We 
take the chessboard (/3 = 1) as an example to study its band structure (shown in Fig. 4(b)). 
From Fig. 4(b) one sees clearly that some degeneracies (including the usual degeneracies 
at the edge of zone B and the accidental degeneracies for points at F — X) are lifted for 
both the i?-polarization and the iJ-polarization. An absolute band gap appears at where 
the degeneracies are lifted at the edge points of zone B. To understand this situation, after 
the unit cell of the photonic crystal in Fig. 3(b) changes to a larger one (denoted by the 
dashed line in Fig. 3(b)), the corresponding first Brillouin of this new photonic crystal is 
denoted by Zone C in Fig. 3(c). The transform matrix between zone C and zone A is A^"^ = 



the Appendix. For the point group symmetry of C^y, we can analyze the degeneracy with 
the two mirror symmetry operators ^3 = ai and 0:4 = a2- The eigenstates with the 2 wave 
vectors on the two opposite edges in the Brillouin zone B are 2-fold degenerate states under 
these two symmetry operations. The degeneracy becomes only 2-fold at points X and M 
now (as compared to the 4-fold degeneracy in the band structure associated with the super 
cell) since |A^^| = 2 here. Therefore, degeneracies must disappear at X and M points and 
each group of 4 bands in Fig. 4(a) break at points X and M to form 2 groups with 2 
bands in each group (see Fig. 4(b)). An absolute band gap Au/uc = 0.070 appears at the 
mid-frequency (of the band gap) u = Uc = 0.605(27rc/a). The reason why they have large 
absolute band gaps can be explained by the present theory of super cells. 

Fig. 4(c) is the corresponding gap map when (3 increases from to 1. The situation is 
similar to the case considered in []T^ (the only difference is that here we use square dielectric 




which corresponds to M = det(A^"^) = 2,{G} = {0,b2} = {0,^} according to 
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FIG. 5: The case of the chess board structure, (a) The super cell including 4 unit cells, (b) Th 
symmetry of the super cell is broken when two square rods increase in size and the other two squar 
rods decrease in size, (c) The photonic crystal when the size of the smaller square rods becomes 
zero. The structure becomes the simple square lattice of square rods. 

rods instead of round air holes). To make a map for the actual procedure of the degeneracy 
breaking, we take P as the parameter. The photonic crystal with additional smaller squares 
(0 < /3 < L we call it case 2) is considered in [15] to be the resulted photonic crystal by 
adding the smaller ones to corners of the square lattice of the square rods [jJ = 1, we call it 
case 3). Here, we take both case 2 and case 3 to be the results of changing the translation 
symmetry of the photonic crystal when the additional square is of equal sized (/? = 0, we 
call it case 1). They have the same degeneracy breaking properties with the chessboard 
structure as mentioned above. Thus larger absolute band gaps can be expected similarly by 
choosing an appropriate value of l3. When P > 0.76, an absolute band gap appears around 



□ 



15 



I I TM-mode 
I I TE-mode 
H Both 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



P 



FIG. 6: The gap map for Fig. 5 with 0</3<l(l — /3is the side length ratio of the smaher 
squares to the larger ones). 

ijj = 0.6(27rc/a). It is more useful to use Au/uc to describe the PBGs due to the scaling 
property of a photonic crystal. In the band structure associated with the super cell, Auj/uc 
remains almost unchanged. A maximum Au/uc = 0.071 occurs when f3 = 0.93. 

Following the same procedure, from a chessboard photonic crystal shown in Fig. 5(a) we 
can obtain a photonic crystal formed with square dielectric rods of two different sizes (see 
Fig. 5(b)) . The ratio of the side lengths between the smaller rods and the larger rods is 
1 — (3 with < (3 < 1. Fig. 6 gives the gap map when the filling ratio is fixed to / = 0.35 
with < f3 < 1 and the inclusion material has a dielectric constant e = 11.4. From Fig. 
6, one sees that there is no absolute band gap for the chessboard case (when /3 = 0) in 
the frequency range of < u < 2iTc/a. When (3 = 1, the structure becomes the simple 
square lattice of square rods of the same size (see Fig. 5(c)), which has an absolute band 
gap Auj/uJc = 0.0453 with the mid-frequency uJc = 0.7231(27rc/a). It is thus not surprising 
that with appropriate parameters an absolute photonic band gap exists for the 2D square 



lattice of square dielectric rods as considered in |jTT|. The maximal gap of Alu/lUc = 0.0717 
(much larger than the case of inclusions with single size) occurs a.t u = Uc = 0.7449(27rc/a) 
when the ratio of the size lengths for the two inclusion rods is (3 = 0.31. 

A triangular lattice of air columns (see Fig. 7(b)) has been found to have a large absolute 
band gap. The rule of thumb can be employed to give a reasonable explanation. Here we 
explain how a gap appears from the view of changing the translation symmetry. Although 
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a 

< ► 




FIG. 7: The case for a rectangular lattice of air holes. The ratio of the two side lengths is 
(a) The super cell including 4 unit cells, (b) The symmetry of the super cell is broken when 
two diagonal air holes reduce size to zero and the other two diagonal air holes increase in size. 
The symmetry- broken super cell is marked by the thick solid lines, (c) The first Brillouin zone A 
(marked by the dashed lines) associated with the original unit cell and the first Brillouin zone B 
(marked by the solid lines) associated with the super cell of Fig. 7(a). The first Brillouin zone C 
associated with the new unit cell (marked by the dashed hnes of Fig. 7(b)) of the new photonic 
crystal of Fig. 7(b) is marked by the dotted lines. F, X, J, M are the symmetry points. 

the triangular lattice has a high symmetry, it can be viewed as the result of symmetry 
breaking from a super cell of another photonic crystal shown in Fig. 7(a). Prom Figs. 8(a) 
and 8(b) we can see clearly how the degeneracies are lifted at the edge points X and M (cf. 
Fig. 7(c)) and an absolute band gap is created in the band structure associated with the 
super cell when the symmetry of the super cell is broken. 
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FIG. 8: Band structures calculated with the plane wave expansion method (with 961 plane waves). 
The relative permittivity for the background medium is e = 13. The solid lines denote the H- 
polarization case and the dashed lines are for the E-polarization case, (a) The band structure 
associated with the super cell of Fig. 7(a). The radium of the air holes is r = 0.5a. (b) The 
band structure associated with the symmetry-broken super cell (marked by the thick solid line 
in Fig. 7(b)) The filling factor is / = 0.836. The absolute band gap l^uj/uc = 0.169 appears at 
mid-frequency lOc = 0.4936(27rc/a). 

As a final numerical example, we break the symmetry of a super cell shown in Fig. 9(a) 
by changing both the sizes (for all the nine square rods) and the positions (except the central 
square rod) of the square rods, however, with the dielectric filling factor / fixed. Fig. 9(b) 
is a resulted structure. The band structure with parameters pi — 0.194a, p2 = 0.236a, p3 = 
0.374, / = 0.5132 are shown in Fig. 10. From this figure one sees that the degeneracies 




18 




FIG. 9: Symmetry breaking of a super cell by changing both the sizes and positions of the inclusions, 
(a) The super cell including 4 unit cells, (b) The symmetry of the super cell is broken when both 

the sizes (for all the 9 square rods) and the positions (except the central square rod) of the square 
rods are changed (however, with the dielectric filling factor / is fixed). 

are lifted at the edge points and two absolute band gaps are created at higher normalized 
frequencies, namely, a large gap AuJi = 0.072(27rc/a) at cUc = 1. 291 {'In c/ a) and another gap 
Au}2 = 0.043 at ujc = 1.142(27rc/o). Note that it is easier to fabricate a photonic crystal 
with the absolute band gap occurring at a higher normalized frequency. 

IV. CONCLUSION 

In the present paper, we have presented a new way of explaining or creating photonic 
band gaps through analyzing the degeneracy of the band structure associated with a super 
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FIG. 10: The band structure associated with the symmetry-broken super cell shown in Fig. 9(b) 
with / = 0.5132 and e = 11.4. The solid lines denote the H-polarization case and the dashed lines 
are for the E-polarization case. The other parameters are pi = 0.194a,p2 = 0.236a, = 0.374a. 
Through the symmetry-breaking of the super cell, one obtains a band structure with two absolute 
band gaps, namely, one gap Aa;i = 0.0724 at ui = 1.2905(27rc/a), and another gap Alo2 = 0.0427 
at W2 = 1.1424(27rc/a). 

cell. The band structure associated with a super cell of a photonic crystal has degeneracies 
at the edge of the first Brillouin zone if the photonic crystal has some kind of point group 
symmetry. We have analyzed these degeneracies and presented two theorems on degeneracies 
in the band structure associated with the super cell. These theorems and analysis are valid 
in any dimensional space (and for any polarization in the 2D case) and does not require the 
investigation of the field distribution. Photonic band gaps can be created through lifting 
these degeneracies by changing the transform matrix between the super cell and the unit cell 
which is caused by changing the translation group symmetry of the photonic crystal. Many 
numerical examples have been given in the present paper to illustrate this. In the 2D case, 
the El-polarization and the H-polarization have the same properties of degeneracies. The 
existence of photonic band gaps for many known 2D photonic crystals have been explained 
through the degeneracy analysis of the band structure associated with the super cell. Some 
photonic crystal structures with large or multiple band gaps have also been found by breaking 
the symmetry of the super cell. 
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Appendix: The proof of the Lemma and a method to find the set {G} 



We can take 3 kinds of primary transformations for the integer matrix while keeping 
the absolute value of det(A^) unchanged. The first transformation is to multiply a column 
or row by ±1. The second transformation is to interchange two columns or rows. The third 
transformation is to add one column or row with k times of another column or row {k is an 
integer). Each primary transformation can be expressed with a left or right multiplication of 
the matrix by an integer matrix. Furthermore, the determinants of these integer matrices 
(associated with these primary transformations) are 1 and their inverses are also integer 
matrices. Under these primary transformations the determinent of the matrix is kept 
unchanged and the resulted matrix A^' is still an integer matrix. 

Here we give a specific procedure for obtaining a diagonal matrix N' = D hj taking 
these transformations. First, we interchange the column or the row with the second kind of 
transformation to make A^^n the minimal among A'^ij and A^ji , i = 1, 2, 3. Then we make Nn 
positive (if it is negative) with the first kind of transformation. If all the integers and 
N^-j^ {i = 1, 2, 3) are divisible by A''{^, we can make all and N--^ (except N[^) zero by using 
the third kind of transformation. If any of N[- or N--^ is not divisible by N[-^, its column or 
row can be subtracted with k times of A^^n so that the remaining component at or N-^ 
position becomes smaller than N[i. The column or row is then interchanged with the first 
column or the first row (with Nn) to make smaller. This can be done repeatedly until 
all and A^^^'^ are divisible by A''^^^ or A''(;^ = 1. Then we can make all and A^^^'^ but A''^^^ 
to be zero. Applying a similar way to N22, we can make A'23 = = 0- Therefore, the 
integer matrix A^^ is diagonalized as 



where D,P are all integer matrices. Here P ^ is the inverse of P and det(P ^) = 
det{Q) = 1. Thus, we have det(L') = ^16^2^3 = det(A^^) = M. 
From Eq. (5) we have 



N'^ = P~^DQ, 
Id, \ 



(Al) 



D 



^2 



V rfs y 



3 3 



(A2) 



i=i j=i 
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_ _ _ 3 _ 3 

We can define two other basic vectors b and b' by bj = J2 Pij^j ^-^id b • = X) Qij^'j- Below 

i=i i=i 
we will find {G} in terms of the vectors b^. It follows from Eq. (14) that 

h,=j2D^,^=j:6,,d,^. (A3) 

j=l 3=1 

Since 

3 3 3 3 

^n^bj = riiP^^ijhj = J^^'j^J ^ {^nihi}, (A4) 

i=l *ii=l j=l i=l 

3 3 3 3 

i=l i,j=l j=l «=1 

3 _ 3 

we see that the set {J2 ^ibi} = {Z^ riibi} {ni,n\ are arbitrary integers). Similarly, we can 

i=l i=l 

show that 

{X:nM} = {X:n,b',}. (A6) 

i=i 1=1 

3 _ 3 

Therefore, the sets {G} and {G'} can be written as{G|G = rijbj = njdjb^} and 

{G'|G'= En,b;}. 

1=1 

Since any integer ni can be written as = lidi + rrii with < mj < — 1, we see that the 
_ _ 3 _ 

set {G|G = J2 n^i^i, n^i — 0, di — 1}. Obviously, we have {G} n {G} = and there are 

i=l 

did2ds — M possible combinations of (mi, m2, ms) for G (i.e., the set {G} has M elements). 

3 _ _ 3 _ 

Also, for any two of them, e.g. Gi = ^Mi, G2 = J2 ^'i^'ij < rrii < di — 1,0 < m'^ < 

i=l i=l 

__3_3_3 

di — 1, the difference Gi — G2 = iriih'^ — Y "m't^i — Y {"nii — "m'd^'i does not belong to 

i=l 1=1 i=l 

{G} since < jm^ — < di — 1. 

Obviously, any G' e {G'} can be written as 

3 3 
G' = 5:n,b: = ^(W,b: + m,b:), 

1=1 1=1 

= G + G, (A7) 

where G G {G} and G G {G}. 

As a numerical example in the 2D case, we can consider the super cell indicated by 
the thicken line in Fig. 3(b) where the unit cell is indicated by the dashed line (the su- 
per cell includes 2 unit cells). Obviously, we have a'l = ai + a2 and a'2 = ai — a2. 
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Then we have N'^ — 




Following the above procedure, we obtain N'^ — 



. D 




1 
1 1 



Q 




. Therefore, we have di = 1,^2 = 2, and {G} = {0,^} = 

{0, b2}. Since Gi — G2 = h'2 in this special example and we have the following general form 
for G - nibi+n2b2 = ni(b;+b'2)+n2(-b;+b^) = (ni-n2)b; + (ni+n2)b^. G1-G2 e {G} 
will end up with ni = 77,2 = 1/2 which is contradictory to the requirement that ni and n2 
are integers. Therefore, for this special example we also see that Gi — G2 ^ {G}. 
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